Enhanced excitonic effects in the energy loss spectra of LiF and Ar at large 
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It is demonstrated that the bootstrap kernel [3| for finite values of q crucially depends upon the 
matrix character of the kernel and gives results of the same good quality as in the q ^ limit. 
The bootstrap kernel is further used to study the electron loss as well as absorption spectra for Si, 
LiF and Ar for various values of q. The results show that the excitonic effects in LiF and Ar are 
enhanced for values of q away from the F-point. The reason for this enhancement is the interaction 
between the exciton and high energy inter-band electron-hole transitions. This fact is validated by 
calculating the absorption spectra under the influence of an external electric field. The electron 
energy loss spectra is shown to change dramatically as a function of q. 
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One of the promising routes for tailoring electro- 
magnetic interactions in systems is to make use of 
excitons- excitonic condcnsation[^ cxciton-plasmon 
interaction [3, [Hj and strong coupling of excitons with 
other inter-band optical transitions. A lot of applications 
of excitonic manipulation have found their way into sur- 
face physics, nano-structures and nanotubcs. Extended 
solids, however, have not been treated as serious candi- 
dates for such effects. This is mainly due to the fact 
that excitonic effects are in general (but not exclusively) 
stronger in lower dimensional systems @. 

As for solids, much of the attention of excitonic 
physics has been focused on the determination of opti- 
cal absorption spectra (specially in the long wavelength 
{q 0) limit). Determination of optical absorption 
spectra in this limit requires an accurate description of 
electron-hole effects, which in turn requires a compu- 
tationally expensive many-body treatment at the level 
of the Bethe-Salpeter equation (BSE)0-[ll- Alterna- 
tively, this electron-hole physics can also be effectively 
treated using time-dependent density functional theory 
(TDDFT)[l5[. However, there exist only a few exchange- 
correlation (xc) kernels within TDDFT which are capa- 
ble of accurately describing the excitonic effects arising 
from strong electron-hole interactions; the computation- 
ally demanding but accurate nano quanta kernel 3, 12 1 
which is derived from the BSE, the long-range corrected 
(LRC) kernel[ll,[il[ll] which has the form a/|gp (with 



a being a system dependent external parameter) and the 
recently proposed bootstrap kernel [ij. The main feature 
which makes these kernels capable of accurately treat- 
ing the absorption spectra (in g — > limit) is their 
l/|gp dependence H^, a feature necessary to capture the 
electron- hole physics. 

In contrast to this the optical absorption spectra and 
the electron energy loss spectra (EELS) at finite values 
of q (away from the F-point) , are known to be accurately 



treated |2lU22| hy the adiabatic local density approxima- 
tion (ALDA) j23j. which does not have the depen- 
dence. This then raises an interesting question about the 
validity of the kernels which accurately treat the q — s- 0, 
for finite values of q. 

In the present work, the EELS and absorption spec- 
tra for Si, diamond, LiF and Ar are calculated, us- 
ing the bootstrap kernel [2| implemented within the ELK 
codefi^- This choice is motivated by the fact that the 
bootstrap kernel is computationally not demanding and 
it does not depend upon any system-dependent param- 
eters. For prototypical materials with bound exlectron- 
hole pair, LiF and Ar, there exists a strong coupluing 
between exciton with high energy inter-band transitions, 
which leads to strong excitonic effects away from the F- 
point. This coupling provides a useful handle to manip- 
ulate the excitonic physics, since inter-band transitions 
can easily be modified by changing external parameters 
like pressure and/or electric field. 

The bootstrap xc-kcrnel within TDDFT for calculating 
linear response reads: 



-i(q,w = O)w(q) e-^{c[,uj = 0) 
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where, v is the bare Coulomb potential, xo is the response 
function of the non-interacting Kohn-Sham system and 
eo(q, i^) = 1 — i;(q)xo(q7w) is the dielectric function in 
the random phase approximation (RPA) . All these quan- 
tities arc matrices in the basis of reciprocal lattice vectors 
G. With a view towards studying the excitonic proper- 
ties of solids for any value of q, we first show that the 
bootstrap kernel for finite values of q leads to accurate re- 
sults for the absorption spectra as well as the EELS. This 
is demonstrated in Fig. [1] which compiles results for the 
dielectric response and EELS of a small bandgap insula- 
tor, Si, a medium badgap insulator, diamond, and a large 
bandgap insulator, LiF. The choice of these materials is 
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motivated by the abundance of available experimental 
data and because they have been used as prototypical 
test cases for the influence of many-body corrections. 
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FIG. 1: Upper two panels- imaginary part of the dielectric 
tensor as a function of energy in eV for Si. Results obtained 
using the bootstrap kernel, RPA and ALDA are compared 
with experimental data from Ref. [2^ . The top panel con- 
tains results for |q| — 0.53a.u. and the middle panel for 
|q| — 0.80a. It. parallel to the [111] direction. The bootstrap 
results are obtained in two different ways; (red line) the kernel 
is a full matrix in the reciprocal space with |G|max = 4a. ?i. 
and (pink dotted line) only the head of the kernel is used 
i.e. |G|max = Oa.it. The bottom panel contains results for the 
EELS, for different values of q (indicated in the figure) as a 
function of energy in eV for LiF and diamond. The results 
obtained using the bootstrap kernel are shown with full line, 
the experimental data (from Ref. [25h with dots and the BSE 
results (also taken from Ref [25l ) with dashed line. The re- 
sults for different values of q are shifted vertically for clarity. 

The upper two panels of Fig. [1] are presented the re- 
sults for £(q, w) of Si. The results obtained using the 
bootstrap kernel are compared to the results obtained us- 
ing the ALDA and to the experimental data. It is clear 
that as far as the absorption spectrum is concerned the 
bootstrap kernel (treated as a full matrix in reciprocal 
space) gives results in good agreement with experiments. 



Interestingly, for small frequencies the results obtained 
using the bootstrap kernel give an upper bound to the 
experimental data, while the results obtained using the 
ALDA give a lower bound. The RPA results (also shown 
in Fig. [T|), as expected, totally miss the excitonic physics, 
which in this case shows up as the shifting of the spectral 
weight to lower frequencies. 

In the q — !• limit the G = G' = component of 
fxc is the most important one and hence the bootstrap 
procedure can be thought of as a self-consistent method 
for obtaining the system dependent parameter a of the 
LRC. While, for finite values of q the matrix character of 
fxc is crucial and the bootstrap kernel is significantly dif- 
ferent from the LRC kernel. The importance of including 
higher G vectors in i^c is demonstrated in Fig. [TJ results 
for Si show that the dielectric function obtained using 
only the G = G' = component of are much higher 
in magnitude and in relatively poor agreement with the 
experimental data. 

EELS corresponds to the negative of the imaginary 
part of eQQ{q,uj), where 00 stands for the G = G' = 
component of the dielectric tensor. The lower panel of 
Fig. [T]contains results for the EELS of LiF and diamond. 
Three different values of q in the T — X direction are 
presented for LiF. Within the first BZ the experimental 
data and BSE results show three main peaks which are 
well reproduced by the bootstrap kernel. On going from 
0.23 to 0.48 rX the plasmonic peak at 25 eV gets smaller 
in magnitude, a feature which is again well captured by 
the bootstrap kernel. Experiments, BSE and bootstrap 
results show that outside the first BZ (for q = 1.50TX) 
EELS is highly suppressed. These results indicate that 
the bootstrap kernel captures the change in — Im[e~^] as 
a function of q very well. We note, however, that the 
magnitude of the peaks is slightly overestimated by the 
bootstrap kernel and, for small energies, the peaks are 
blue shifted by ^1 eV compared to experiment. This 
shifting of the excitonic peaks to higher frequencies and 
the overestimation of their magnitude was also a feature 
of the absorption spectra in the long wavelength limit. 

For diamond, the magnitude of the EELS obtained us- 
ing the bootstrap kernel is overestimated compared to the 
experiment. A similar overestimation is also seen in the 
BSE resultst25j|. In fact, we find that the results obtained 
using the bootstrap kernel are in very good agreement 
with the BSE results. As in the case of Si, for LiF and 
diamond as well, the RPA results are shifted to higher 
frequencies and are missing excitonic physics. 

It is apparent from the examples above that the boot- 
strap results for values of q away from the F-point, have 
the same good agreement with experiments as those in 
the q — >■ limit. With this in hand we can now use the 
bootstrap kernel to study the excitonic effects in various 
directions in the BZ. 

In Fig. [2] the results for e(q, w) are shown for LiF 
and Ar for various values of q. The upper panel con- 
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FIG. 2: Dielectric tensor obtained using the bootstrap kernel 
as a function of energy in eV. Upper panel contains results for 
LiF and the Lower panel for Ar in the T — L — X direction. 
The upper-left most panel also contains results at the F-point 
for LiF exposed to an external electric field. 



tains results for LiF in the T — L — X direction[26j. It 
is immediately clear from this that the excitonic peak in 
LiF becomes stronger (spectral weight moves to lower fre- 
quency) away from the F-point; on going from F to L, the 
excitonic peak (^13-14 eV) becomes larger in magnitude 
and at the same time the inter-band transition (IT) peak 
^20 eV diminishes and moves to lower frequencies (from 
20 eV to 18 eV). In the L — X direction the excitonic and 
IT peaks remain almost the same. 

Similar results are also seen for Ar; on going from F 
to L, the spectral weight moves towards the excitonic 
peak (~12eV) which becomes stronger and shifts to lower 
frequencies. This is accompanied by a steady diminishing 
of the IT peak around 18 eV. Beyond 0.6FL. both the 
excitonic and the IT peaks are unchanged. The RPA 
for both LiF and Ar misses this excitonic physics for all 
values of q. 

The above results point towards a strong coupling be- 
tween the excitonic and the IT peaks; as the excitonic 
peak gains weight, the IT peak diminishes. This obser- 



vation could be used in the future to tune excitonic effects 
via manipulation of inter-band transitions. In order to 
validate this hypothesis we have performed calculations 
of the q — >■ absorption spectrum for LiF in presence 
of an external electric field. This field is an artifice used 
to lift the degeneracy and split the bands. This split in 
the bands causes the IT peak at 20 eV to diminish and 
move to lower energies. The excitonic peak gains weight 
(spectral weight moves to lower energy) and the results 
for e in the q — >■ limit closely resemble those of the L- 
point (see top-first panel of Fig. [5]). Such a tuning of the 
excitons has been performed before: a coupling of exci- 
tons and surface plasmons was used to enhance excitonic 
effects in low dimensional systems [3, [13 ■ 
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FIG. 3: Electron energy loss spectrum (— /m[eQQ^(q, cj)])in 
arbitrary units as a function of energy in eV for LiF. Results 
obtained using the bootstrap kernel are shown with full lines 
and the results obtained using the RPA are shown with dashed 
lines. Upper panel shows results in the F — L direction and 
lower panel in the X — L direction. Results for various values 
of q are shifted vertically for clarity. 



eQQ'(q, w) is an important quantity, not just because 
it can be directly compared to experiments, but also be- 
cause it is required for an accurate determination of the 
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screened Coulomb interaction, W = e~^v. In Hcdin's 
theoretical foundation of many-body perturbation the- 
ory the effective interaction, W, is a crucial concept. 
It is needed as an essential ingredient if one wishes to 
sum Feynman diagrams to obtain the Green's function 
of a system0, 0. In Fig. [3] we present results for the 
EELS of LiF as a function of q and frequency. The EELS 
changes dramatically as a function of q; the peak at 16 
eV gains height on moving in the T — L direction (up- 
per panel), while at the same time the peak at 25 eV 
diminishes in magnitude and broadens. This leads to the 
EELS for q — ^ being significantly different from that at 
q = [0.5, 0.5, 0.5]. A similar sharpening of the peak at 16 
eV is also seen in the L — X direction (lower panel); A 
broad double peak at -^18 cV changes to one sharp peak 
at 16 eV and a small shoulder at ~20 eV. In the L — X di- 
rection (unlike the T — L) the peak at 25 eV stays almost 
unchanged as a function of q. Similar strong changes in 
EELS as a function of q, in the T — U direction for LiF, 
have been reported before in Ref. . 

The RPA results (also shown in Fig. [3]) not only miss 
the low energy excitonic effects |28|. but also show very 
little variation as a function of q in both the T ~ L and 
X — L directions. Most many-body calculations use RPA 
dielectric functions for screening the Coulomb potential. 
However, it is clear from the present study that RPA di- 
electric functions can wildly differfrom those which prop- 
erly include electron- hole physics. TDDFT with the 
bootstrap kernel is a computationally efficient method 
to accurately determine £~^(q, w) for use in many-body 
perturbation theory [l^. 

To summarize, in the present work it is demonstrated 
that results obtained using the bootstrap kernel for finite 
values of q have the same accuracy as those in the q — 5- 
limit. We use this kernel to make the predictions that in 
the prototypical materials LiF and Ar excitonic effects 
arc enhanced away from F. The reason for this enhance- 
ment is attributed to the interaction between the exci- 
ton and other high energy inter-band transitions. This 
observation is reinforced by the application of an exter- 
nal electric field and noting an inverse proportionality 
in strengths between the inter-band and excitonic peaks. 
It is further demonstrated that the EELS also changes 
dramatically as a function of q. This strong change in 
EELS is missing within the RPA, and hence it is highly 
desirable to use a TDDFT dielectric function to screen 
the Coulomb interaction as input for many-body pertur- 
bation theory. 
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